arXiv:1509.00004v2 [astro-ph.GA] 24 Sep 2015 


Draft version September 25, 2015 

Preprint typeset using style emulateapj v. 5/2/11 


THE r-PROCESS IN MAGNETOROTATIONAL SUPERNOVAE 
Takuji Tsujimoto^ and Nobuya Nishimura^ 

^National Astronomical Observatory of Japan, Mitaka-shi, Tokyo 181-8588, Japan; taku.tsujimoto@nao.ac.jp 
^Astrophysics Group, Keele University, ST5 5BG Keele, UK 
Draft version September 25, 2015 

ABSTRACT 

One of the hottest open issues involving the evolution of r-process elements is fast enrichment in 
the early Universe. Clear evidence for the chemical enrichment of r-process elements is seen in the 
stellar abundances of extremely metal poor stars in the Galactic halo. However, small-mass galaxies 
are the ideal testbed for studying the evolutionary features of r-process enrichment given the potential 
rarity of production events yielding heavy r-process elements. Their occurrences become countable 
and thus an enrichment path due to each event can be found in the stellar abundances. We examine 
the chemical feature of Eu abundance at an early stage of [Fe/H] < — 2 in the Draco and Sculptor 
dwarf spheroidal (dSph) galaxies. Accordingly, we constrain the properties of the Eu production in 
the early dSphs. We find that the Draco dSph experienced a few Eu production events, whereas Eu 
enrichment took place more continuously in the Sculptor dSph due to its larger stellar mass. The 
event rate of Eu production is estimated to be about one per 100-200 core-collapse supernovae, and 
a Eu mass of ~ (1 — 2) x 10~^Mq per single event is deduced by associating this frequency with the 
observed plateau value of [Eu/H] ^ —1.3 for [Fe/H] > —2. The observed plateau implies that early Eu 
enrichment ceases at [Fe/H] « — 2. Such a selective operation only in low-metallicity stars supports 
magnetorotational supernovae, which require very fast rotation, as the site of early Eu production. We 
show that the Eu yields deduced from chemical evolution agree well with the nucleosynthesis results 
from corresponding supernovae models. 

Subject headings: nuclear reactions, nucleosynthesis, abundances — galaxies: dwarf — galaxies: evo¬ 
lution — galaxies: individual (Draco, Sculptor) — stars: abundances 


1. INTRODUGTION 

What is the astrophysical object producing r-process 
elements and how has r-process enrichment proceeded 
through cosmic time? Clear answers to these questions 
are still veiled. Recent updates from both observational 
and theoretical studies of the r-process site, in particu¬ 
lar by the detection of a near-infrared light bump in the 
afterglow of a sh(3 rt-duration y-ray burst (sGRB)(e.g., 
Tanvir et al.||2013 ) and successful nucleosynthesis calcu¬ 


lations (e.g., Rosswog et al. 2014 Wanajo et al.| 2014 


Goriely et al.l^Oibp , have indicated neutron star (Nb) 


mergers as the most promising astronomical site. Core¬ 
collapse supernovae (CC-SNe) driven by neutrino heating 
are another candidate for the r-process site, buy they are 
unable to synthesize heavy r-process elements such as Eu 
and Ba. The required conditions for r-process nucleosyn¬ 
thesis are not achieved in the proto-NS wi nds due to the 
strong effect of neutrino absorption (e.g. , Thompson et 
al.|2001 Fischer et al.|2010| Wanajo|2013 ). 'This tact has 
turned out to strengthen the INS merger origin scenario. 

New research on r-process synthesis has renewed inter¬ 
est in the Galactic chemical evolution (GCE) of r-process 
elements in terms of the enrichment taking place through 
NS mergers. Though previous researchers have pointed 
out the fata l problems with the notion of enrichment by 
NS mergers dMathews fc Cowan 1990 Argast et al. 2004), 
Tsujimoto &: Shigeyama 12014) successfully reproduced 
an observed Eu feature iii the Galaxy using a GCE model 
that considers a unique propagation of the ejecta of NS 
mergers in the scheme of hierarchical galaxy formation. 
They also found that an early release of Eu from (at least 


some fraction of) NS mergers with a short timescale of 
^ 10^ year is an inevitable condition for explaining the 
presence of very low-metallicity ([Fe/H] —3) stars en¬ 

riched by r-process elements. Such a prompt enrichment 


1 Matteucci et al. 

2014 

Ishimaru et al. 

2015 

Wehmeyer 

et al. 

2015), and 

can be a possible channel 

in popula- 

tion s 

ynthesis models ( 

Belczynski et al. 

2006 

Dominik 


merger depends critically on the binary evolution during 
a common envelope phase which is poorly understood. 
In addition, the above requirement is in contrast with 
a typical NS merger time of 3-4 Gyr deduced from the 
analysis of Swift sGRB samples (Wanderman & Piran 


2015^ 


i'sujimoto & Shigeyama (2014) have extended the dis- 
cussion of the site of r-process elements to nearby dwarf 
spheroidal (dSph) galaxies. They have found compelling 
evidence for NS mergers as the origin of r-process el¬ 
ements, in the observed chemical feature of Eu in the 
Draco, Sculptor, and Carina dSphs; these dSphs ex¬ 
hibit a constant [Eu/H] of ~ —1.3 with no apparent 
increase in the Eu abundance over a metallicity range 
of — 2 [Fe/H] — 1. This implies no Eu production 
events as more than 10^ CC-SNe increase the galactic Fe 
abundance. This level of rarity of Eu production is com¬ 
patible with the frequency e xpected for NS mergers (Tsu¬ 
jimoto & Shigeyama 2014). The next question is how 
these dSphs which do not undergo any Eu production 
events for [Fe/H] ^ —2 are enriched up to [Eu/H]^ —1.3 
until [Fe/H]~ —2. Tsujimoto et al. (2015) examined Eu 
abundances for [Fe/H]< —2 in the Draco and Sculptor 
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Fig. 1. — Observed Eu feature for the Draco (red circles), Sculp¬ 
tor (blue circles), and Carina (crosses) dSphs. A brief illustration 
of early enrichment in the Draco and Sculptor dS hs is shown. The 
observ ed data are updated from those of Fig. 3 in |Tsujimoto et al.| 
l|2015|l including [Eu/H] deduced from [Ba/H] for stars with upper 
limits on Eu data or no detection of Eu together with new data. 
The Eu floor is assigned at [Eu/H]~ —5, which is implied by th e 


ll2ute , iKirb v &: CohenI (|2012jl7Tsimon et ah ) 2U15t' andij a blonk a 
et al. fcuiht lor the 'Scul ptor dSph: and Ish etrone et al.l ||2UU3||. 
Gemasle et al.| l|2012[l, and [Venn et ai.|||201^[| tor tne Uarina dapA. 


dSphs from a purely observational viewpoint and con¬ 
cluded that [Eu/H] remarkably increases from < —4 to 
^ —1.3 inside these dSphs, including a conspicuous in¬ 
termittent increase in the Draco dSph. These revealed 
Eu features suggest that the contributor of Eu in the 
early dSphs is active only in the low-metallicity environ¬ 
ment with an occurrence rate that is not as frequent as 
in CC-SNe. 

Besides the mechanism of CC-SNe by neutrino heat¬ 
ing in the standard scenario, another proposed option 
for r-process enrichment is an explosion t riggered by fast 
rotat ions and high magnetic helds (e.g., Takiwaki et al. 
20091. These magnetorotational CC-SNe (MR-SJNe) can 
produce heavy r-process elements and thus may con- 
tribute to their enrichment during early star formation 


(|Fuiimot oetalj2008| Wi nteler et al. | [20T2 Nishimura et 


al.||2015 ' VVehmeyer et al.||2015 ). in addition, it is possi- 

ble to presume that the einergence of MR-SNe is inclined 
toward very low-metallicity stars in which the rotationa l 
velocity is expected to be high (e.g., Meynet et al.|2006 ). 

In this Letter, we provide compelling evideiice tor 
heavy r-process enrichment by MR-SNe at an early stage 
of star formation inside dSphs. This Letter is organized 
as follows. Theoretical arguments on the chemical fea¬ 
ture of Eu abundance in the Draco and Sculptor dSphs 
are first presented in Section 2. Then, the interpreta¬ 
tion of these arguments is incorporated into the models 
of chemical evolution, and the results are shown in Sec¬ 
tion 3. Subsequently we discuss the astrophysical site of 
early Eu production in terms of nucleosynthesis in MR- 
SNe (Section 4) and finally conclude with Section 5. 


2. EARLY EU ENRICHMENT IN THE DRACO AND 
SCULPTOR dSph GALAXIES 

Chemical features of Eu abundance for faint (i.e., less 
massive) dSphs in the luminosity range of 3 x 10® — 2 x 
10®L©, i.e., the Draco, Carina, and Sculptor dSphs, is di¬ 


vided into two discrete trends: (i) a remarkably increas¬ 
ing [Eu/H] for [Fe/H] — 2, and (ii) a broadly constant 
[Eu/H] for [Fe/H] — 2. The first feature is also sup- 
ported by the trend of [Ba/H] increasing versus [Fe/H] 
( Tsujimoto et al.|[2015 ). It suggests that Eu production 
episodes occur only in the low-metallicity regime. Here 
we attempt to unravel how the early Eu enrichment pro¬ 
ceeds in the Draco and Sculptor dSphs (note that the 
Carina dSph has only one datum for [Fe/H]< —2). The 
observed correlation of [Eu/H] with [Fe/H] is shown in 
Figure 1. In this figure, for the stars ([Fe/H]< —2) mea¬ 
sured with only upper limits of Eu abundance or with¬ 
out the detection of Eu, we deduce [Eu/H] from their 
Ba abundances with an assumed pu re r-process ratio of 
[Ba/Eu]=—0.89 (Burris et al. 2000). This procedure is 
validated by the fact that the stars with [Fe/H] ^ —2 
exhibit a [Ba/Eu] ratio compatible with the value as¬ 
sumed in the translation (i.e., [Ba/Eu]=—0.89) and thus 
are considered to be of a pure r-pro cess origin before 


s-pro cess nucleosynthesis can occur (Tsujimoto et al. 


20151. 

First we see the enrichment path in the Draco. It 
begins with a jump of [Eu/H] from < —4.6 to ~ —2.5 
between metallicities [Fe/H]« —3 and —2.5. To pin 
down the metallicity of the first Eu production episode 
more precisely, acquisition of more [Eu/H] data between 
—3 <[Fe/H]< —2.5 are awaited. This includes verifica¬ 
tion of the possibility of intrinsic scatter in [Eu/H] at 
[Fe/H]f« —3 implied by two data points. Though the 
subsequent path leading to a plateau of [Eu/H]~ —1.3 
is hard to detect, it is possible to predict the presence 
of a few jump-like increases in [Eu/H] as indicated by 
dotted arrows if we consider the errors in [Fe/H] and 
[Eu/H] of individual stars. In fact, such a small num¬ 
ber of Eu production events is strongly supported by a 
relatively high metallicity where the first Eu jump hap¬ 
pens. The corresponding metallicity is determined by 
how many CC-SNe producing Fe occur before the first 
event, which is directly connected to the frequency of Eu 
production event. For the Draco dSph with metallici¬ 
ties of [Fe/Hjfti —3 to —2.5, this gives about three to six 
events in total as shown by our calculations in the fol¬ 
lowing section. Here we will deduce the actual rate of 
events from this total number of 3-6. The stellar metal- 


et al. 

2011 

) indicates that ^ 387o of stars populate the 

metal 

iicitv 

ran 

2 :e of 

Fe/HJ< —2. Using the Kroupa IMF 
), we obtain the corresponding stel- 

(Kroupa et al. 

1993 


lar mass for this m etallicity range of ^ 

(Martin et al.|2008l and Hnd that this stellar population 
presumably hosted ~ 600 total CC-SNe in total. Ac¬ 
cordingly, we conclude that r-process production events 
happen at a rate of one per ^ 100 — 200 CC-SNe in the 
early evolution of Draco. In addition, we can estimate 
the average Eu yield of a single event. The average metal¬ 
licity of [Fe/H]^ —1.9 ( Kirby et al.||2011 ) is expected to 
result from the conversion mass fraction of the initial 
gas to stars of ^ 10%, which is calculated by models of 


dwar f galaxies (Tsujimoto|2011 Tsujimoto & Shigeyema 
2012[ see also s3). Then, the hypothesis that an initial 


gas cloud with a mass of 3.2 x 10 ®Mq was enriched to 
[Eu/H]=—1.3 by three to six events leads to a Eu yield 
of ~ (1 - 2) X 1O-®M0. 

























































































3 


In the Sculptor dSph, a rise from the Eu floor hap¬ 
pens at [Fe/H]~ —3.5, a lower metallicity than Draco’s, 
i.e., [Fe/H] —3. Then, this Eu jump is followed by 
a continuous increase in [Eu/H] if we ignore an outlier 
showing a very low Eu abundance ([Eu/H]=—2.75) at 
[Fe/H]~ —2.4. We note that three unpublished d ata sets 
of [Ba/H] at [Fe/H]« —2.4 ( Tolstoy et al.|[200^ suggest 
([Eu/H]) = —2.2. Such a relatively continuous feature is 
naturally understood if we take into consideration that 
the observed stellar mass of Sculptor is about eight times 
more than that of the Draco. We expect that the oc¬ 
currence frequency of early Eu production events in the 
Sculptor will increase by the same factor, and as a result, 
more than 20 occurrences of Eu production are likely to 
build the smooth enrichment path as observed. Irrespec¬ 
tive of its higher event rate, the final enrichment level 
should be the same as for Draco since the mass ratio of 
initial gas between the two dSphs is broadly equivalent 
to the ratio of the Eu event rate between the two. There¬ 
fore, different mass-scaled dSphs, i.e., Draco, Carina, and 
Sculptor, share the same plateau of [Eu/H]~ —1.3. How¬ 
ever, the evolutionary path to the plateau varies in ac¬ 
cordance with the number of Eu production episodes. 

3. CHEMICAL EVOLUTION MODELS 

To validate the theoretical interpretation presented in 
the previous section, we calculate the evolution of Eu 
abundance for the Draco and Sculptor dSphs, adopt¬ 
ing the following properties of Eu producers. We as¬ 
sume two production sites for Eu: NS mergers and one 
that emerges selectively in a low-metallicity regime, likely 
MR-SNe (see Section 4). Due to the se vere rarity of 
NS mergers (one per thousands of CC-SNe,|Tsujimoto &' 
Shigeyani^|2014 )), the two faint dSphs did not host any 
MS mergers over their whole lifetimes. MR-SNe, how¬ 
ever, contribute to Eu enrichment in the early phase with 
a moderately rare event rate. Here we assume that MR- 
SNe operate with a rate of one per 200 CC-SNe in the 
metallicity range of [Fe/H] < —2 for the reference case. 
This rate gives three occurrences in Draco as already 
discussed, while 21 total events for Sculptor are deduced 
from its stellar mass and MDF (Kirby et al. 2011). In our 
calculations, each MR-SN is assumed to occur after ev¬ 
ery 200 CC-SNe with a progenitor metallicity lower than 
[Fe/H]=—2. For the Eu mass ejected from a single MR- 
SN, we adopt the yield of 2 x By contrast, for 

the first MR-SN for both dSphs, a smaller amount of Eu 
mass is found to give a better fit to the observed jumps 
from an Eu floor, and thus 7 x 1O“®M0 is assumed as 
the Eu mass of the first MR-SN. As discussed in §4, both 
Eu yields are within the predictions from hydrodynami- 
cal simulations of MR-SNe. In addition to the reference 
case, for Draco we calculate another case in which the 
MR-SN rate is assumed to be one per 120 CCSNe with 
a Eu yield of 1.2 x \Q~^Mq (4 x \Q~^Mq for the first 
MR-SN). 

The chemical evolution models are modified from the 
GCE model to account f or the star formation rate (SFR ) 
and the IMF as done by|Tsujimoto & Shigeyema|(|2012|). 
To reproduce the observed IVlDh's with a mean [Ee/H] 
of —1.93 and —1.68 for the Dr aco and Sculptor dSphs, 
respectively (Kirby et al. 2011), we adjust the SFR to¬ 
gether with a steep IMF (i=— 1.6) as in the Large Magel¬ 
lanic Cloud. Here we assume that the period of star for- 



Fig. 2. — Predicted evolutions of Eu abundance for the Draco 
(red curves) and the Sculptor (blue curve) dSphs, compared with 
the observed data shown in Fig. 1. The Eu production events are 
assumed to occur at a rate of one per 200 CC-SNe (solid curves) or 
one per 120 (350) CC-SNe (dashed red (blue) curve) for [Fe/H] < 
—2 while staying dormant for [Fe/H] > —2. A resultant smaller 
event rate (three: solid curve; five: dashed curve) in the Draco is 
attributed to its lower stellar mass compared with the Sculptor. 
The pre-enriched level is set at [Eu/H] = — 5. 


mat ion continues for 4 Gyr for Draco (cf. | Dolphin||2002 1 
and 6 Gyr for Sculptor (jde Boer et al.|2012|), respectively. 
The initial Eu abundance (i.e., the Eu floor) is set to be 
[Eu/H]=—5. The results from our models are shown in 
Figure 2. 

The observed values of [Fe/H] corresponding to the 
jump from a pre-enriched level of [Eu/H] for the two 
dSphs are well reproduced by our models. The higher 
jump of [Eu/H] at higher [Fe/H] in Draco can be inter¬ 
preted as a result of a lower frequency of MR-SNe owing 
to a smaller total stellar mass than Sculptor. The subse¬ 
quent ladder-like increasing feature of [Eu/H] in Draco is 
also in good agreement with observations. On the other 
hand, a gradually increasing trend after the first jump 
until the plateau-like abundance (i.e., [Eu/H]^ —1.3) 
is predicted by the model for Sculptor. Three data 
points giving the low upper limits of [Eu/H] between 
—3.5 <[Fe/H]< —3.3 suggest a more delayed first jump. 
Accordingly, the result of another model in which 12 to¬ 
tal events are assumed to occur with a rate of one per 350 
CCSNe is shown by the dashed blue curve. The agree¬ 
ment of our models with observations for both Draco and 
Sculptor suggest that Eu production episodes occur more 
frequently than NS mergers, but much less frequently 
than CC-SNe and with a high Eu yield. 


4. NUCLEOSYNTHESIS IN MR-SNE 

As discussed in the previous sections, the early chemi¬ 
cal evolution of the Draco and Sculptor dSphs demands 
a large amount of Eu production such as ~ (1 — 2) x 
1Q~^Mq per single event. This value far exceeds the 
assumed amount (10“^ — 1Q~^Mq) from the ejecta of 
canonical C C-SNe with neutrin o-driven winds in GCE 
models (e.g., Argast et al.|2004 |, while it is smaller than 
the expectation from an NS merger (~ 1O“^M0). Here 
we claim that the implied Eu yields agree well with 
the predictio n from MR-SNe, b ased o n the recent cal¬ 
culations by Nishimura et al. (2015). In their mod¬ 
els, a typical case (i.e., Bll/tl.OO) results in a Eu mass 
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Fig. 3.— Nucleosynthesis abundance patterns predicted by MR- 
SN models for three different cases, co mpared with the solar abun¬ 
dance pattern (|Arlandini et al.||1999|). See the text for the nucle¬ 
osynthesis models. 


M(Eu) = 0.724 X 1O“^M0 in the jet-like ejecta and a 
total mass of r-process nuclei of 1.53 x Inside 

the ejecta, the values of Yg vary between 0.15 and 0.5 and 
an average in the region where r-process operates (i.e., 
Ye < 0.3) is (Ye) = 0.192. The resultant nucleosynthesis 
pattern is shown by the red curve in Figure 3, compared 


Arlandini et al. 

1999 

tion rates basec 

on t 


1999). In this calculation, we used reac- 


which is available in the reactio n rate library REACLIB 
(Rauscher & Thielemann 2000), leading to a relatively 
low Ba/Eu ra tio of |Ba/Eu|= — 1.44. However, a recent 
improvement ( Moller el al.||2012 ) indicates that the pro¬ 
duction of isotopes in the region incl uding the rare eart h 
peak has increased by a few factors (Kratz et al. 2014), 
which may lead to a match with the solar pattern, giv¬ 
ing M(Eu) = 1.63 X 1Q~^Mq a nd a [Ba/Eu1 rati o close 
to [Ba/Eu]=—0.89. In addition, Nishimura et al. (2015) 
pointed out that the total mass of ejecta can be larger by 
a few factors in more energetic jet-like explosion models. 
This strong jet model gives [Ba/Eu] = —0.55 (the blue 
curve). Accordingly, the Eu mass of ~ (1 — 2) x 
is indeed within the prediction by MR-SN models. 

Besides the above jet-like explosions, another scenario 
that would eject a significant amount of Eu is from 
MR-SNe, which release a massive amount of moderate 
neutron-rich ejecta. According to recent updated magne¬ 
tohydrodynamic si mulations that resolve the magnetoro- 
tational instability (Sawai & Yamada 2014), the ejecta 
due to neutrino heating enhanced by fast rotation and 
large magnetic fields have much lower values of Yj, for 
the cases of explosion without jet-like outflows. As a re¬ 
sult, the ejecta are a moderate neutron-rich state. To 


obtain the nucleosynthesis results for the corresponding 
case, we set higher Y^ values inside the ejecta in the jet 
model (see Section 5.2 of Nishimura et al. 2012) and 
obtain results for (Yg) = 0.254, as shown by the dashed 
curve in Figure 3. The ejecta mass i ncluding the to¬ 


tal r- process nuclei of 2.15 x 10 '^Mq (Nishimura et al. 
2015) gives M(Eu) = 0.85 x which is broadly 


equivalent to the mass predicted by the jet-like explosion 
models. 


5. CONCLUSIONS 

We show that the properties of r-process production 
events that are important for early star formation, i.e., 
the frequency and yield, can be assessed from the re¬ 
cently revealed Eu feature of two less massive dSphs. 
The frequency is estimated to be about one per 100-200 
CC-SNe in a low-metallicity regime, while this frequency 
decreases to a level similar to or less than that of NS 
mergers for [Fe/H] > —2. We identify high frequency 
events at low-metallicity together with the implied high 
Eu yield as MR-SNe. Accordingly, we have the following 
picture of r-process enrichment among galaxies. In dwarf 
galaxies such as the Draco and Sculptor dSphs, MR-SNe 
at an early phase are the only contributors of the nu¬ 
cleosynthesis of heavy r-process elements. However, in 
more massive galaxies, including our own, in addition 
to MR-SNe with their metallicity-dependent frequencies, 
NS mergers enrich the ISM with heavy r-process elements 
over the whole lifetime of a galaxy at a rate of one per 
thousands of CC-SNe. Since enrichment by MR-SNe is 
expected to start much earlier than NS mergers, the pres¬ 
ence of extremely metal-poor stars enriched by r-process 
elements as observed in the Galactic halo is naturally 
explained by our scheme. 

We predict that the frequency of MR-SNe is fairly 
small in the relatively metal-rich galaxies. Thus their 
occurrences should be detected with a high rarity in the 
present-day Universe. MR-SNe are likely to be associ¬ 
ated with magnetars, the rate of w hich is estimated t o 
be 0.003 — 0.06 yr“^ in the Galaxy (Muno et al. 2008). 
We accordingly anticipate that a few percent or less of 
magnetars are the end result of MR-SNe. A link between 
MR-SNe and rare ph enomena such as super luminous SNe 
j Dess art et al.|2012[) and ultra-long GRBs (Greiner et al. 
2015| ), which may Ije energized by magnetars, is worth 
investigating in future work. 
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